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Results from a phenomenological analysis of W and Z hard diffractive hadroproduction at high 
energies are reported. Using the Regge factorization approach, we consider the recent diffractive 
parton density functions extracted by the HI Collaboration at DESY-HERA. tn addition, we take 
into account multiple Pomeron exchange corrections considering a gap survival probability factor. 
It is find that the ratio of diffractive to non-diffractive boson production is in good agreement with 
the CDF and DO data. We make predictions which could be compared to future measurements at 
the LHC. 



PACS numbers: 12.40.Nn, 13.85.Ni, 13.85.Qk, 13.87.Ce 



I. INTRODUCTION 

Recently, the diffractive processes are attracting much 
attention as a way of amplifying the physics programme 
at proton colliders, including new channels searching for 
New Physics. The investigation of these reactions at high 
energies gives important information on the structure of 
hadrons and their interaction mechanisms. Hard diffrac- 
tive processes, such as the diffractive production of mas- 
sive electroweak bosons and dijets, allow the study of the 
interplay of small- and large-distance dynamics within 
QCD. The existence of a hard scale provides the normal- 
ization of the Born term diagram. For boson hadropro- 
duction, single-diffractive dissociation can occur charac- 
terized by the existence of one large rapidity gap, which 
can be represented by the Pomeron exchange. At high 
energies, there are important contributions from unita- 
rization effects and the suppression of the single-Pomeron 
Born cross section due to the multi-Pomeron contribu- 
tions depends, in general, on the particular hard process. 
At the Tevatron energy, yfs = 1.8 TeV, the suppression 
is in the range 0.05-0.2 [H, 0,1,13, whereas for LHC 
energy, y/s — 14 TeV, the suppression appears to be of 
order 0.08-0.1 [3, [2, El] • Therefore, the correct treatment 
of the multiple scattering effects is crucial for the relia- 
bility of the theoretical predictions of the cross sections 
for these diffractive processes. 

In the present study, our motivation to perform a new 
analysis on diffractive boson production is twofold: pro- 
duce updated theoretical estimations compatible with 
the current Tevatron data on single diffractive W and 
Z hadroproduction 0, Q and to perform reliable predic- 
tions to the future measurements at the LHC. In order 
to do so, we use Regge factorization (single Pomeron ex- 
change) and the corresponding corrections for multiple- 
Pomeron scatterings. Factorization for diffractive hard 
scattering is equivalent to the hard-scattering aspects of 
the Ingelman and Schlein model pj], where diffractive 
scattering is attributed to the exchange of a Pomeron, 
i.e. a colorless object with vacuum quantum numbers. 
The Pomeron is treated like a real particle, and one con- 



siders that a diffractive electron-proton collision is due 
to an electron-Pomeron collision and that a diffractive 
proton-proton collision is due to a proton-Pomeron col- 
lision. Therefore, the diffractive hard cross sections are 
obtained as a product of a hard-scattering coefficient, a 
known Pomeron-proton coupling, and parton densities in 
the Pomeron. The parton densities in the Pomeron have 
been sistematicaly extracted from diffractive DIS mea- 
surements. In particular, the quark singlet and gluon 
content of the Pomeron is obtained from the diffractive 
struture function F®^ 1 (sjp, P, Q 2 )- Recently, a new anal- 
ysis of these diffractive parton distributions has been pre- 
sented by the HI Collaboration in DESY-HERA. On 
the other hand, it is well known that the single Pomeron 
approach produces results that overestimate the exper- 
imental values by a large factor [t| [nj. Thus, in the 
present analysis the corresponding multiple Pomeron ex- 
change corrections will be taken into account. 

This paper is organized as follows. In the next sec- 
tion, we present the main formulae to compute the inclu- 
sive and diffractive cross sections for W and Z hadropro- 
duction. We show the details concerning the parame- 
terization for the diffractive partons distribution in the 
Pomeron, extracted recently in DESY-HERA. In addi- 
tion, we present the theoretical estimations for the gap 
survival propability factor that will be used in the com- 
parison of our results with experimental measurements 
from Tevatron and extrapolations to the LHC energy. In 
the last section we present our numerical results and per- 
form predictions to future measurements in CERN LHC 
experiment. The compatibility with data is analysed, 
possible additional corrections are investigated and the 
comparison with other approaches are considered. 



II. DIFFRACTIVE HADROPRODUCTION OF 
MASSIVE GAUGE BOSONS 

Let us start by introducing the main expressions 
to compute the inclusive and diffractive cross sections. 
For the hard diffractive processes we will consider the 
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Ingelman-Schlein (IS) picture 0], where the Pomeron 
structure (quark and gluon content) is probed. The start- 
ing point is the generic cross section for a process in which 
partons of two hadrons, A and B, interact to produce a 
massive electroweak boson, A + B — > (W ± /Z°) + X, 



da 



dx a dxb 



^2 fa/ A (Xa , fJ. 2 ) fb/B ( x b , 



da(ab -> [W/Z] X) 
di 



plicit formulation is described in terms of Regge theory. 
Therefore, we can rewrite the parton distribution as 



Xafa/A{Xa, M ) = / dx W / d/3 



''J dtf w/A (x w , t) 

X Pfa/AP, ), (1) 

\ XjpJ 



where Xifi/h(xi, jj, 2 ) is the parton distribution function of 
a parton of flavour i — a, b in the hadron h = A, B. The 
quantity da/dt gives the elementary hard cross section 
of the corresponding subprocess and p 2 = is the 

hard scale in which the pdf 's are evolved in the QCD evo- 
lution. Equation above express the usual leading-order 
QCD procedure to obtain the non-diffractive cross sec- 
tion. Next-to-leading-order contributions are not essen- 
tial for the present purposes. 

In order to obtain the corresponding expression for 
diffractive processes, one assumes that one of the 
hadrons, say hadron A, emits a Pomeron whose par- 
tons interact with partons of the hadron B. Thus the 
parton distribution x a f a /A(x a , V?) in Eq. (pj is replaced 
by the convolution between a distribution of partons in 
the Pomeron, (3f a /w(Pj M 2 ), an d the "emission rate" of 
Pomcrons by the hadron, fw/h(xTP,t). The last quan- 
tity, fw/h{xw,t), is the Pomeron flux factor and its ex- 



and, now defining f(xjp) 
tains 



I°oo dt f tp /k{xw,t), one ob- 



Xafa/A^Xa, ^) 



/ dx w f(xjp) — f a /T?{ — ,M 2 )- ( 2 ) 
J x-sp Xtp 



Concerning the W diffractive production, one consid- 
ers the reaction p+p(p) — > p+ W(— > e v)+ X, assuming 
that a Pomeron emitted by a proton in the positive z di- 
rection interacts with a p (or a p) producing that 



subsequently decays into 



The detection of this re- 



action is triggered by the lepton (e + or e~) that appears 
boosted towards negative rj (rapidity) in coincidence with 
a rapidity gap in the right hemisphere. 

By using the same concept of the convoluted structure 
function, the diffractive (single diffraction, SD) cross sec- 
tion for the inclusive lepton production for this process 
becomes 



da SB 

uu lepton _ 



drj e 



a . b 



dxw 
x w 



f(xjp) / dE T f a/w (x a , n )fb/p(p)(x b , M ) 



V a \ G% 



6 s fi 



(3) 



where 



and 



Xb 



M w 



M w e-v- 



A±y/(A*-1) 



a t V(a 2 -i) 



t = -E T M w A + y/(A 2 - 1) 



p? = M w . The quantity V a b is equal to the Cabibbo- 
Kobayashi-Maskawa matrix element if e a + e/, = ±1 and 

(4) zero otherwise, where a, b denote quark flavors and e q the 
fractional charge of quark q. The upper signs in Eqs. (j4|) 
and © refer to W + production (that is, e + detection). 
The corresponding cross section for W~ is obtained by 

(5) using the lower signs and t <-> u. 

In a similar way, the cross section for the diffractive 
hadroproduction of neutral weak vector boson Z is given 
by 

(6) 



with A = Mw/2E~t, Et being the lepton transverse 
energy, Gp is the Fermi constant and the hard scale 



J 



4 D (v^) = E /— /— / — /(*-)/./-(*.. 

. J X W J Xb J X a 



P 2 ) fb/p(p)(x b , /i 2 ) 



^C z ah G F M 2 z 



da(ab ZX) 
di 



(7) 
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where = 1/2 - 2\e q \ sin 2 9 W + A\e q \ 2 sin 4 9 W , with 
9w being the Weinberg or weak-mixing angle. The def- 
initions for x a ,b are similar as for the W case and now 
/j 2 = M§. The values of the electroweak parameters that 
appear in the various formulae were taken from the par- 
ticle data handbook [ll[, and we use only four flavors 
(u, d, s, c) in the weak mixing matrix, with the Cabibbo 
angle 9 C = 0.2269. 



A. The Pomeron Flux Factor 

An important element in the calculation of hard 
diffractive cross sections is the Pomeron flux factor, in- 
troduced in Eq. (fT|). We take the experimental analysis 
of the diffractive structure function Q , where the Xjp de- 
pendence is parameterised using a flux factor motivated 
by Regge theory [l2[, 



/n>/p(3:ip>*) = a tp 



2a P (t)-l 



(8) 



where the Pomeron trajectory is assumed to be linear, 
ajp(t) = Q!np(0) + a.' w t, and the parameters B-p and a' w 
and their uncertainties are obtained from fits to HI FPS 
data fl3j . The normalisation parameter is chosen 
such that Xtp ■ j t mm /ip/p dt = 1 at = 0.003, where 
|i min | ~ m 2 p x\ j (1 — x^p) is the minimum kinematically 
accessible value of |f|, m p is the proton mass and |f C ut| 
1 .1) GeV 2 is the limit of the measurement. 

The flux factor above corresponds to the standard 
Pomeron flux from Regge phenomenology, based on the 
Donnachie-Landshoff model [l4|. On the other hand, 
there is an alternative Pomeron flux, proposed first by 
Goulianos [ljl, which considers it as a probability den- 
sity. Thus, the integral over the diffractive phase space 
could not exceed the unit and the standard flux should 
be normalized. For instance, see Refs. [To| for previ- 
ous phenomenology using the normalized flux in boson 
hadroproduction. 



B. The Pomeron Structure Function 

In the estimates for the diffractive cross sections, we 
will consider the diffractive pdf's recently obtained by 
the HI Collaboration at DESY-HERA [8]. The Pomeron 
structure function has been modeled in terms of a light 
flavour singlet distribution E(z), consisting of u, d and 
s quarks and anti-quarks with u = d=s = u = d = s, 
and a gluon distribution g(z). Here, z is the longitudi- 
nal momentum fraction of the parton entering the hard 
sub-process with respect to the diffractive exchange, such 
that z = [3 for the lowest order quark-parton model pro- 
cess, whereas < (3 < z for higher order processes. The 
quark singlet and gluon distributions are parameterised 




FIG. 1: (Color online) The rapidity distribution of elec- 
tron and positron generated in inclusive and diffractive W 
hadroproduction at ^/i =1.8 TeV (see text). 



at Qq with the general form, 



zfiiz, Qq) 



(1 - z) Cl exp 



0.01 



L (1-*) 



(9) 



where the last exponential factor ensures that the diffrac- 
tive pdf's vanish at z = 1. The charm and beauty quarks 
are treated as massive, appearing via boson gluon fusion- 
type processes up to order a 2 . To determine experimen- 
tally the diffractive pdf's, the following cuts have been 
considered: j3 < 0.8, M x > 2 GeV and Q 2 < 8.5 GeV 2 , 
mostly in order to avoid regions influenced by higher 
twist contributions or large theoretical uncertainties [f| . 

For the quark singlet distribution, the data require 
the inclusion of all three parameters A q , B q and C q in 
equation [5] By comparison, the gluon density is weakly 
constrained by the data, which are found to be insen- 
sitive to the B g parameter. The gluon density is thus 
parameterised at Qq using only the A g and C g param- 
eters. With this parameterisation, one has the value 
Ql = 1.75 GeV 2 and it is referred to as the 'HI 2006 
DPDF Fit A'. It is verified that the fit procedure is not 
sensitive to the gluon pdf and a new adjust was done 
with C g = 0. Thus, the gluon density is then a simple 
constant at the starting scale for evolution, which was 
chosen to be Q§ = 2.5 GeV 2 and it is referred to as the 
'HI 2006 DPDF Fit B'. The quark singlet distribution is 
well constrained, with an uncertainty of typically 5 — 10% 
and good agreement between the results of both fits [8| . 



C. The Gap Survival Factor 

In the following analysis we will consider the sup- 
pression of the hard diffractive cross section by multi- 
Pomeron scattering effects. This is taken into account 
through a gap survival probability factor. There has 
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been large interest in the probability of rapidity gaps 
in high energy interactions to survive as they may be 
populated by secondary particles generated by rescatter- 
ing processes. This effect can be described in terms of 
screening or absorptive corrections, which can be esti- 
mated using the quantity [T(ij |: 



TABLE I: Data versus model predictions for diffractive W 
hadroproduction (cuts -ET min = 20 GeV and i r < 0.1). 



<|S| 2 >= 



J\A(s,b)\ 2 e- n( - s ^ d 2 b 
J\A(s,b)\ 2 d 2 b 



(10) 



where A is the amplitude, in the impact parameter space, 
of the particular process of interest at center-of-mass en- 
ergy y/s. The quantity is the opacity (or optical den- 
sity) of the interaction of the incoming hadrons. This 
suppression factor of a hard process accompanied by a 
rapidity gap depends not only on the probability of the 
initial state survive, but is sensitive to the spatial distri- 
bution of partons inside the incoming hadrons, and thus 
on the dynamics of the whole diffractive part of the scat- 
tering matrix. 

For our purpose, we consider two theoretical estimates 
for the suppression factor. The first one is the work of 
Ref. 0] (labeled KMR), which considers a two-channel 
eikonal model that embodies pion-loop insertions in the 
pomeron trajectory, diffractive dissociation and rescat- 
tering effects. The survival probability is computed for 
single, central and double diffractive processes at sev- 
eral energies, assuming that the spatial distribution in 
impact parameter space is driven by the slope B of the 
pomeron-proton vertex. We will consider the results for 
single diffractive processes with 2B = 5.5 GeV -2 (slope 
of the electromagnetic proton form factor) and with- 
out N* excitation, which is relevant to a forward pro- 
ton spectrometer (FPS) measurement. Thus, we have 

< \S\ 2 >kmr= 0.15 for y/s = 1.8 TeV (Tevatron) and 

< |S| 2 >kmr= 0.09 for y/s = 14 TeV (LHC). 

The second theoretical estimate for the gap factor is 
from Ref. [TtJ (labeled GLM), which considers a sin- 
gle channel eikonal approach. We take the case where 
the soft input is obtained directly from the measured 
values of er tot , a e i and hard radius Rh- Then, one has 

< \S\ 2 >glm= 0.126 for y/s = 1.8 TeV (Tevatron) and 

< |S1 2 >glm= 0.081 for y/s = 14 TeV (LHC). We quote 
Ref. [TtJ for a detailed comparison between the two ap- 
proaches and further discussions on model dependence of 
inputs and consideration of multi-channel calculations, ft 
should be stressed that our particular choice by KMR and 
GLM (single channel) models is in order to indicate the 
uncertainty (model dependence) of the soft interaction 
effects. It is worth to mention that some implementa- 
tions of GLM model include the results of a two or three 
channel calculation for < IS"! 2 >, which are considerably 
smaller than the one channel result [l3l ■ 



III. RESULTS AND DISCUSSION 

In the following, we present our predictions for hard 
diffractive production of Ws and Z's based on the pre- 





Rapidity 


Data (%) 


Estimate (%) 


1.8 TeV 


M < 1.1 


1.15±0.55 [5] 


0.715 4= 0.045 


1.8 TeV 


r? e | < 1.1 


1.08 4=0.25 [6] 


0.715 4=0.045 


1.8 TeV 


1.5 < |r? e | < 2.5 


0.64 4= 0.24 [6] 


1.7 4= 0.875 


1.8 TeV 


Total W -> ev 


0.89 4= 0.25 [6] 


0.735 ± 0.055 


14 TeV 


M < 1 




31.1 ± 1.6 



vious discussion. These predictions are compared with 
experimental data from Refs. [H, [|| in Tables I and II. 
In addition, estimations for the LHC are presented. In 
the numerical calculations, we have used the new HI pa- 
rameterizations for the diffractive pdf's Q. The 'HI 2006 
DPDF Fit A' was considered and one verifies that the re- 
sults are not quite sensitive to a replacement by 'HI 2006 
DPDF Fit B'. For the usual pdf's in the proton (anti- 
proton) we have considered the updated MRST2004F4 
parameterization [l8l |. which is a four-fixed-flavour ver- 
sion of the standard MRST2004 parton distributions. As 
the larger uncertainty comes from the gap survival factor, 
the error in the predictions correspond to the theoretical 
band for < \S\ 2 >. In the theoretical expressions of pre- 
vious section only the interaction of pomerons (emitted 
by protons) with antiprotons (protons in LHC case) are 
computed, that means events with rapidity gaps on the 
side from which antiprotons come from. The experimen- 
tal rate is for both sides, that is events with a rapidity 
gap on the proton or antiproton side. Therefore, we have 
multiplied the theoretical prediction by a factor 2 in order 
to compare it with data. 

Let us start by the diffractive W production. In order 
to illustrate our investigation, in Fig. [T] we present the 
rapidity distribution of the electron (dot-dashed lines) 
and positron (solid lines) generated in both inclusive and 
diffractive hadroproduction in Tevatron for y/s = 1.8 
TeV. The diffractive cross sections are not corrected by 
gap survival factor and they are given byEq. ©. In this 
case, the diffractive production rate is approximately 7 
% (using the cut \r/\ < 1) being very large compared 
to the Tevatron data. When considering the gap survival 
probability correction, the values are in better agreement 
with data. When considering central W boson fraction, 
-1.1 < Ti e < 1.1 (cuts of CDF and DO @, Q), we ob- 
tain a diffractive rate of 0.67 % using the KMR estimate 
for < | SI 2 >, whereas it reaches 0.76% for the GLM esti- 
mate. The average rate considering the theoretical band 
for the gap factor is then R w = 0.715 ± 0.045 %. This 
result is consistent with the experimental central values 
i?cp F = 1.15 % and i?™ = 1.08 %. The agreement 
would be better whether the sub-leading reggeon con- 
tribution is added, which was not considered in present 
calculation. In Ref. [ijjj], it was shown that its intro- 
duction considerably enhances the diffractive ratio in the 
Tevatron regime. Considering the forward W fraction, 
1.5 < |?7 e | < 2.5 (DO cut), one obtains R w = 0.83 % 
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TABLE II: Data versus model predictions for diffractive Z { 



hadroproduction (cuts -ET mil 


, = 16 GeV and 


if < 0.1). 


t/s Rapidity 


Data (%) 


Estimate (%) 


1.8 TeV Total Z -> e+e" 


1.44 ± 0.80 [6] 


0.71 ±0.05 


14 TeV Total Z e+e" 




30.26 ± 1.41 



for KMR and R w = 2.58 % for GLM, with an averaged 
value of Rw = 1-7 ± 0.875 %. In this case, our estimate 
is larger than the central experimental value R™ — 0.64 
%. For the total W -> ei/ we have i? w = 0.68 % for 
KMR and R w = 0.79 % for GLM and the mean value 
Rw = 0.735 ± 0.055 %, which is in agreement with data 
and consistent with a large forward contribution. Fi- 
nally, we estimate the diffractive ratio for the LHC en- 
ergy, y/s — 14 TeV. In this case we extrapolate the pdf's 
in proton and diffractive pdf's in Pomeron to that kimc- 
natical region. This procedure introduces somewhat ad- 
ditional uncertainties in the theoretical predictions. We 
take the conservative cuts |77 e | < 1, -Er mJn = 20 GeV for 
the detected lepton and x ff < 0.1. We find Rw — 32.7 % 
for KMR gap survival probability factor and Rw = 29.5 
% for GLM, with a mean value of i?^ 10 = 31.1 ± 1.6 
%. This means that the diffractive contribution reaches 
one third, or even more, of the inclusive hadroproduc- 
tion even when multi-Pomeron scattering corrections are 
taken into account. The reason of this enhancement is 
the increasingly large diffractive cross section. The re- 
sults presented above are summarized in Table I. The 
experimental errors have been summed into quadrature. 

Now, we present the investigations for the diffractive 
Z hadroproduction. When the gap survival factor is not 
considered the diffractive cross section is given by Eq. 
([7]), producing a diffractive rate of 6.2 %. This value is 
once again higher than the Tevatron data by a factor five. 
When considering the gap survival correction, we verify 
an agreement with experiment. For the total Z — > e + e~ 
we obtain a diffractive rate of 0.66 % using the KMR 
estimate for < \S\ 2 >, whereas it reaches 0.76 % for the 
GLM estimate. The average value gives Rz = 0.71±0.05, 
which is consistent with the experimental result R™ = 
1.44 + 0.61 — 0.52. A rough extrapolation to LHC energy 
gives R z = 31.67 with KMR gap factor and R z = 28.85 
for GLM, with a mean value #| HC = 30.26 ± 1.41. We 
again consider the conservative cuts £r min = 16 GeV and 
ip < 0.1. This estimate follows similar trend as for the 
W case. The results presented above are summarized in 
Table II. The experimental errors have been summed into 
quadrature. 

Our results can be compared with previous calcula- 
tions in diffractive boson hadroproduction. For instance, 
in Refs. [To| one uses IS approach with a normalized 
Pomeron flux [l5| and the corresponding diffractive pdf's. 
The data description for the W case is reasonable. How- 



ever, the calculations are only compared to the CDF [5j 
data and they are somewhat larger than ours. In Ref. 
[l9| a hard Pomeron flux is considered, i.e. a ff (0) ~ 1.4, 
and multiple scatterings are taken into account by a 
Monte Carlo calculation. In addition, for Tevatron ener- 
gies the reggeon contribution is added. The results are 
compared only to CDF data [5] for the W production 
and the description is consistent with experiment. It is 
interesting the fact that a hard Pomeron flux could mimic 
the multi-Pomeron suppression or the effect of normal- 
izing the standard Pomeron flux. Finally, we need call 
attention to the uncertainty in the determination of the 
gap survival probability. The estimates considered here 
(KMR and GLM) are compatible with each other for 
the case of single diffractive processes. However, recent 
calculations using one channel eikonal model give larger 
values for < |Sf > @,Hj. For instance, in Ref. [H an 
eikonal QCD model with a dynamical gluon mass (DGM) 
was considered. Using a gluon mass m g — 400 MeV, one 
obtains < |S| 2 >dgm (Tevatron) = 27.6 ± 7.8 % and 
< \S\ 2 >dgm (LHC) = 18.2 ± 7.0 %. These values give 
R w (y^ = 1.8 TeV) ~ 1.23 % and R z {*/s = 1.8 TeV) ~ 
1.21 %. This illustrates the size of uncertainty when con- 
sidering different estimates for the gap probability. 

In summary, we have shown that it is possible to ob- 
tain a reasonable overall description of hard diffractive 
hadroproduction of massive gauge bosons by the model 
based on Regge factorization supplemented by gap sur- 
vival factor. For the Pomeron model, we take the recent 
HI diffractive parton density functions extracted from 
their measurement of F% . The results are directly de- 
pendent on the quark singlet distribution in the Pomeron. 
We did not observe large discrepancy in using the dif- 
ferent fit procedure for diffractive pdf's (fit A and B). 
We estimate the multiple interaction corrections taking 
the theoretical prediction of distinct multi-channel mod- 
els, where the gap factor decreases on energy. That is, 
<|S'| 2 >~ 15 — 17.5 % for Tevatron energies going down 
to < |S| 2 >~ 8.1 - 9 % at LHC energy. We find that the 
ratio of diffractive to non-diffractive boson production is 
in good agreement with the CDF and DO data when con- 
sidering these corrections. The overall diffractive ratio 
for y/s — 1.8 TeV (Tevatron) is of order 1 %. In addi- 
tion, we make predictions which could be compared to 
future measurements at LHC. The estimates give large 
rates of diffractive events, reaching values higher than 
30 % of the inclusive cross section. 
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